Achieving steady-state entanglement of remote micromechanical oscillators by 

cascaded cavity coupling 
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In this paper, we propose a scheme for generating steady-state entanglement of remote microme- 
chanical oscillators in unidirectionally-coupled cavities. For the system of two mechanical oscillators, 
we show that when two cavity modes in each cavity are driven at red- and blue-detuned sidebands, 
respectively, a stationary two-mode squeezed vacuum state of the two mechanical oscillators can 
be generated with the help of the cavity dissipation. The degree of squeezing is controllable by 
adjusting the relative strength of the pump lasers. Our calculations also show that the achieved 
mechanical entanglement is robust against thermal fluctuations of phononic environments. For the 
case of multiple mechanical oscillators, we find that the steady-state genuine multipartite entangle- 
ment can also be built up among the remote mechanical oscillators by the cavity dissipation. The 
present scheme does not require nonclassical light input or conditional quantum measurements, and 
it can be realized with current experimental technology. 

PACS numbers: 



I. INTRODUCTION 

Besides fundamental research interests in quantum 
physics [U , realizing quantum effects of macroscopic ob- 
jects is crucial for potential applications in ultrahigh pre- 
cision measurements and quantum information process- 
ing [2 4]. Thanks to the recent achievements in ground- 
state cooling of micromechanical oscillators via optome- 
chanical coupling [5|-|8j], the emerging field of cavity op- 
tomechanics as an interface between mechanical systems 
and optical field has become a unique platform to study 
uantum behavior of macroscopic mechanical systems [3- 
Using well-established quantum optical techniques, 
optomechanics holds the promise to effectively prepare 
and manipulate nonclassical mechanical states. 

Several schemes have been proposed to establish en- 
tanglement between a mechanical element and the driven 
cavity field II 511 or between vibrating membranes or end 
mirrors [16l4l9| | by optomechanics. Apart from short- 
distance mechanical entanglement, remote entanglement 
between two micromechanical oscillators in separated 
cavities can also be entangled via injecting squeezed light 
or conditional quantum measurements [2fj|422| . It was 
also showed that weak mechanical entanglement between 
two distant optomechanical oscillators can be possibly 
achieved merely by optomechanical coupling 23]. The 
entanglement of remote mechanical elements is of impor- 
tance for constructing long-distance quantum communi- 
cation networks [24 |. 

On the other hand, generating quantum states by 
quantum-reservoir engineering has attracted a lot of at- 
tention recently. In this approach, the interaction be- 
tween system and environment is engineered in such a 
way that the system relaxes into a desired state. The re- 
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suiting quantum states are steady, independent of initial 
conditions, and most importantly robust against incoher- 
ent noise. To date, several schemes have been proposed 
to prepare entan gled states of atomic systems by quan- 
tum dissipation [25l43l| and the dissipative creation of 
steady-state entanglement between two separated cold 
atomic ensembles has been experimentally realized [32j . 

In this paper, we consider the generation of steady- 
state entanglement of remote micromechanical oscilla- 
tors (membranes) by cavity dissipation. We at first in- 
vestigate the entanglement between two micromechan- 
ical membrane oscillators in a cascaded cavity system. 
In each cavity, a membrane oscillator is coupled to two 
nondegenerate cavity modes via parametric and beam- 
splitter-like interactions by driving the relevant cavity 
modes on blue- and red-detuned sidebands respectively. 
For negligible mechanical damping, we find that the cav- 
ity dissipation can pull the two distant mechanical oscil- 
lators into a stationary two-mode squeezed vacuum. It 
is also shown that the two-mode entanglement is robust 
against thermal fluctuations when one takes into account 
the mechanical damping. We then extend the two-mode 
mechanical model to the case of multiple mechanical os- 
cillators in an array of cascaded cavities. We show that 
in this system genuine multipartite steady-state entan- 
glement can be built up among the remote mechanical 
oscillators with the help of the cavity dissipation. 

The reminder of this paper is arranged as follows. In 
Sec. II, the model of two cascaded optomechanical system 
is introduced and the steady-state entanglement between 
the mechanical oscillators is investigated in detail. In 
Sec. Ill, we extend the previous model to the case of multi- 
ple mechanical oscillators in an array of unidirectionally- 
coupled cavities and discuss the generation of multipar- 
tite entanglement among multiple mechanical oscillators. 
At last, we give the conclusion in Sec. IV. 
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II. ENTANGLEMENT OF TWO MECHANICAL 
OSCILLATORS 



A. Model and equations 



As schematically shown in Fig.l, we investigate a sys- 
tem consisting of two identical optical cavities connected 
by unidirectional coupling (33j . In each cavity, two driven 
cavity modes are coupled to a vibrating membrane via 
radiation pressure [35| . The role of the membranes 
could also be played by other mechanical systems such as 
trapped clouds of ultracold atoms (36j . After removing 
the carrier photons with filters, the output quantum fluc- 
tuations from the first cavity are directed to the second 
cavity to drive the corresponding cavity modes. With 
the light fields rotating at their driving frequencies, the 
Hamiltonian of the system is given by 



H/h= ]T [S aj a] aj + S bj b% 



+ U} mj CjCj 



3 = 1,2 
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where cij (aj) and bj (6j) (j = 1, 2) are annihilation (cre- 
ation) operators for the cavity modes and Cj (cj) for the 
mechanical modes of the vibrating membranes in each 
cavity. The cavity-laser detunings 5 Z . = oj z . — v z . (z = 
a, b), with uj z . being the cavity resonant frequencies and 
v z . the corresponding driving frequencies. The mechani- 
cal frequencies of the membranes are denoted by uj mj and 
the optomechanical coupling g Zj — h/rrijUj mj ui Zj / L, 
with L being the cavity length and rrij the effective mass 
of the membranes. The amplitudes of the driving lasers 
\£ Zj | = y/2P Zj k Zj /Hv Zj , where P Zj are the powers of the 
pump lasers and k z - the cavity loss rates of the left cavity 
mirrors. 

We expand the quantum operators as oj = + 5oj , 
where o| are the steady-state classical amplitudes and 
Soj the corresponding quantum fluctuation operators. 
By taking into account cavity losses and mechanical 
damping, the classical amplitudes are obtained as z| = 
£z,/(n Zj + iA Zj ) and cj = T, z 9zA^\ 2 /(uj m . + i-f mj ), 
where A Zj — 8 Zj +2g Zj Re(Cj), n Zj are the cavity loss rates 
from the output mirrors on the right of the cavities, and 
r y mj are the mechanical damping rates. Note that here we 
have assumed the cavity loss rates n Zj ^> k Zj such that 
losses from the left cavity mirrors can be neglected. For 
intense driving fields we have |o*| 2 ^> (So^Soj) and the 
Hamiltonian ([1} can be linearized. Then, by dropping the 
symbol "<5" in the fluctuation operators for simplicity of 
notation, the resulting Langevin equations of motion for 
the quantum fluctuations of the cavity and mechanical 



FIG. 1: (a) Schematic setup of two cascaded cavity- 
optomechanical systems. In each cavity, two cavity modes 
are driven by red- and blue-detuned lasers, respectively, and 
the output quantum fluctuations from the first cavity are di- 
rected to the second cavity to drive the corresponding cavity 
modes, (b) Frequencies of the pumps and cavity modes, and 
the dashed arrows represent the unidirectional coupling be- 
tween the cavity modes. 



modes are obtained as 

aj = - (n aj + iA aj )aj - ig aj (cj + cj) + ^2n aj af(t), 
bj = - (k 6j + iA bj )bj ~ ig hj (cj + cj) + ^2l^bf{t), 
Cj = - (7™, + iu mj )cj - ig aj (aj + aj) - ig bj {bj + &j) 

+ y/2^C? (t), (2) 

where the effective optomechanical coupling g Zj — 
\zj\g aj (z = a, b). The noise operators a 1 1 n (i) and 
b™(t) describe vacuum inputs to the first cavity and sat- 
isfy nonzero correlations (a 1 1 n (t)a I 1 n ^(f')) = S(t — t') and 
(6i n (i)6 i 1 nt (t')> = 8{t - t'). The input noise of the second 
cavity, characterized by the operators a™^) and b^it), 
are from the output fluctuations of the first cavity and 
transmission losses in the coupling. When the output 
quantum field of the cavity mode a\ (6i) is used to drive 
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the cavity mode <X2 (62), then one has 

+ Vi^hK n (t), (3a) 

+ v^T^)^(t), (3b) 

where t] z e [0, 1] (z — a, b) accounts for the imperfect 
couplings between the two cavities. The operators a 2 n (i) 
and b™{t) denote the local vacuum noise input to the 
second cavity. The parameter r\ z = 1 corresponds to 
a lossless unidirectional coupling between the two cav- 
ities, whereas r\ z = describes two independent cavi- 
ties. Note here that the exponential factors in the above 
equations result from the differences between the fre- 
quencies of the relevant pump lasers. In addition, cj n (£) 
are noise operators of the mechanical oscillators which 
have nonzero correlations (c™ (£)c™(i')) = n{ h 5(t — t') 

and {cf{t)c l f(t')) = (n 3 th + l)S{t - t'), where the mean 
thermal phonon numbers at temperature T is given by 
fiL — ^ e riuJm j/ kBT — l) -1 , with ks the Boltzmann con- 
stant. 

Now we choose the detunings 

A Ql = -A bl = -u> mi , A a2 = -A b2 = co m2 , (4) 

i.e., the cavity modes a\ and &2 are pumped by lasers 
which are blue-detuned from their resonance frequencies 
by the mechanical frequencies, while the modes 0,2 and 
b\ are driven by pump lasers which are red detuned by 
the same amount, as illustrated in Fig. 1(b). Therefore, 
the pump frequencies v x . should satisfy 

Va-L - v a2 = (w mi +w m2 ), (5a) 

V bl - Vb 2 = -(w mi +UJrn 2 )- (5b) 

With the above choices of detunings, by perform- 
ing the transformations Zj — > Zje~ lA *i l , z™(t) — > 

4 n (i)e _lA ^* (z = a, b), Cj -J- c 3 e'^ m ^\ and cf(t) 
c™(t)e _lCl ' m j*, and neglecting fast oscillating terms pro- 
portional to e ±, ( |J '»i+ u '»!) t ] the Langevin equations ([2]) 
reduce to 

fli = -K ai ai - ig ai c\ + v / 2K ai a 1 i(t), (6a) 
h = -K bl bx - ig bl C! + v / 2K bl b 1 i(t), (6b) 
a-2 = -Ka 2 a 2 - ig a2 C2 - 2^Jr\ a K ai ^a 2 «i + \/2r] a K, a2 d^ (t) 

+ y/2{l -ri a )K aa ai?(t), (6c) 
h = -Kb 2 h - iQb 2 c\ ~ 2y?r) b K bl K b2 bi + ^/2f] b K b2 b^(t) 

+ y/2(l-r, b ) Kba b^(t), (6d) 

ci = -7 mi ci - ig ai a\ - ig bl h + ^/27 mi c' 1 n (t), (6e) 

C2 = -7m 2 c 2 - ig a2 d2 - igb 2 b\ + v / 27 m2 c 2 n (t). (6f) 

It should be noted that for our approximations to be 
valid, we require our system to be in the resolved side- 
band regime, ui mj ^> n Zj , as well as to satisfy uj mj ^> g Zj . 



The above equations show that in each cavity, the me- 
chanical mode is coupled to the cavity modes via effec- 
tive parametric amplification as well as beam-splitter-like 
mixing. While the former interaction leads to photon- 
phonon entanglement and optical amplification, the lat- 
ter is damping the mechanical modes. If the coupling 
strengths satisfy g bl > g ai and g a2 > gb 2 , optical damp- 
ing is dominant over amplification and both mechanical 
oscillators are cooled. 



B. Two- mode mechanical entanglement 

We can equivalently reexpress Eqs. © 
as \ — 2>x + f ln (t), with the vector x — 

(2-ai j Tjai j ^61 j Dbi j %a 2 7 Va 2 > %b 2 ; Vb 2 j %c\ 7 Vc\ > %c 2 7 Vc 2 ) 7 

in terms of the quadrature operators defined as 
x = (o + ot)/V2 and y = -i(o - while / in (t) 

contains the corresponding noise operator contributions. 
The entanglement between the mechanical systems is 
contained in the 12 x 12 correlation matrix a given 
by cry = (xiXj + XjXi)/2- I n steady-steady state, it 
satisfies Za s + <j s Z t = —T>, where D is the noise matrix 
V l0 S(t - t>) = (ftWfit') + ff (0/| n (*))/2. Since 
we are only interested in the entanglement between 
the two mechanical modes, it is enough to consider 
the reduced correlation matrix a 12 related to the two- 
mode mechanical states. It has the simple structure 

°"12 = ( / "s^t a 2 2 ) 7 where <j\ 2 , a\ 2 , and af 2 are 2x2 

matrices containing the autocorrelations of the two 
systems and their cross-correlations respectively. The 
entanglement between the two mechanical modes can 
be quantified with the logarithmic negativity E12 [13], 
which is defined as 

£ 12 =max[0,-hi(2Ci2)], (7) 

where C12 is given in terms or the reduced correlation 
matrix 

C12 = 2- 1 / 2 V / £(a 12 ) - VE(<7i 2 ) - 4det<7i2, (8) 

with £(012) = detfj^ + det(T 2 2 — 2dettT2 2 . 

Solving Eqs.© numerically and using Eq. we can 
investigate the mechanical entanglement in the system. 
Let us first, however, turn to a regime where we can ob- 
tain analytical results. To this end, we consider the cav- 
ity dissipation rates k z . = k, the perfect cavity couplings 
r) z = 1, and the effective optomechanical couplings 

g ai = gb 2 = gi, g a2 = = 52- (9) 

If the cavity dissipation rate is dominating the dynamics 
of the system, i.e. k ^> {gj,J mj n J th }, the cavity modes 
follow changes of the mechanical oscillators adiabatically 
for times t > 1 /ft. In this case we can eliminate the 
cavity modes and find the simple equations of motion for 
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the mechanical modes c, 



Cl(t) = ~{l mi + 7m)ci(t) + y/2 lmi C?(t) + 8*(fi), 

(10a) 

C 2 (t) = -(7m 2 + 7m)c 2 (*) + x/^W^W + c^(t), 

(10b) 

where j m = {g\ — g\)j k is the net optomechanical damp- 
ing rate. The noise operators c™(i) are given by 



c 2 W = — 7=^1 (*) H 7^°i (*), 



(11a) 
(lib) 




12 3 4 
K/271 [ 10 s Hz 



1.00 



and have non- vanishing correlations 



{&f{t)^{t l ))^2^ 1m N m 5{t-t l ) 1 



3 

(S;\t)&f{t')) = 2j m (N m + l)5(t - t'), 
(ci n ^)4 n (t')> = 2j m y/N m {N m + l)6(t - 0, 



(12a) 
(12b) 
(12c) 



with N m = gf /(g% — g\ ) . The above correlations indicate 
that the two mechanical oscillators are effectively cou- 
pled to a broadband quantum reservoir in a two-mode 
squeezed vacuum state [38]. In the absence of the me- 
chanical damping (7 m . = 0), the mechanical oscillators 
will reduce to the state of the reservoir in the long-time 
limit, i.e., the two-mode squeezed vacuum 



|V)i2 = exp(- 



-rc\c\ 



rcic 2 )|0 Cl ,0 c 



(13) 



with the squeezing parameter r = tanh -1 (31/(72) only 
dependent on the relative strengths of the two pump 
lasers. Therefore, the strong mechanical entanglement 
can be built up in principle just by controlling the ra- 
tio of the strengths of the pump lasers. It should be 
pointed that our scheme is quite different from that in 
Ref. [20J which discussed the establishment of the station- 
ary entanglement between two mechanical oscillators by 
injecting externally squeezed light into the cavities. Here, 
instead of creating entanglement in an external source, 
the entanglement between the mechanical oscillator and 
blue-detuned cavity mode is created via the parametric 
interaction in each cavity, and the photon-phonon entan- 
glement is then transferred to the mechanical oscillators 
with the help of the beam-splitter interaction. 

By taking into account mechanical damping, from 
Eq. flTUl) we have the steady-state values (CjCj) — 

(7m"th +7mA f m)/(7m + 7m) an d (ciC 2 ) = J m M m /(jm + 

7 m ), where we have assumed j mi = 7 m2 = 7™, and 
nj h = njr h = nth for simplicity. It is easy to the en- 
tanglement parameter 

Ci2=2 



31(32 - 3i) ~ nj m n t h 



K7„ 



9l 



.9? 



(14) 



FIG. 2: (a) Dependence of the steady-state mechanical entan- 
glement E\2 on the cavity dissipation rate k for different cou- 
pling ratios #2/(71. Other parameters are the mechanical decay 
rate 7 m = 0, the coupling strength 31/2-71 — 0.1 x 10 5 Hz, and 
the unidirectional intercavity coupling efficiency 3 = 0.95. (b) 
The mechanical entanglement as a function of coupling effi- 
ciency 3 for the relative strength (72/31 = 1.2 and the other 
parameters are the same as in (a). 




FIG. 3: Dependence of the mechanical entanglement on the 
coupling 31 for different values of mean thermal phonon num- 
ber n t h, with the cavity decay rate k/2tt = 4 x 10 s Hz, 
the mechanical damping rate 7 m /2-7i = 100 Hz, the coupling 
32 = 1.53i, and unidirectional coupling efficiency 3 = 0.95. 



Clearly, steady-state mechanical entanglement can be 
achieved at non-zero temperature, provided that the 
mean number of thermal phonons satisfy 



- . 9i92 n 3i x 
n t h < (1 )■ 

K7m 32 



(15) 



Given that the couplings g.j are tunable through the 
pump lasers, this condition demonstrates the robustness 
of steady-state entanglement against thermal noise in the 
mechanical systems. 

We next turn to the numerical results from solving 



5 




FIG. 4: Schematic plot of N vibrating membranes trapped in 
cascaded cavities. 



(iL> m ) and identical cavities, we thus find the effective de- 
tunings 



and for the driving frequencies 



-2uj r , 



V\, 2 



2uj r , 



(16) 



(17) 



With the same procedures and approximations as before, 
the Langevin equations of motion for the cavity modes Zj 
and the mechanical modes c, can be obtained and read 



3-1 



K z Zj - ig z c Z j - 2k z y^y^) 



J-s, 



Eqs.©, which allows us to investigate the entanglement 
property in the regime where the adiabatical elimination 
of the cavity modes is invalid. In Fig. 2 the dependence 
of steady-state mechanical entanglement on the cavity 
decay rate n is plotted for different values of 32/31 and 
the mechanical damping 7,,,. = 0. Consistently with our 
analytic results we observe for large cavity decay k» gj , 
the entanglement becomes saturated and independent of 
k. The increase of the entanglement with decreasing cou- 
pling ratios 32/31 is also evident in this regime. Further- 
more, for the converse situation k <C gj, we also observe 
the steady-state entanglement, although to a smaller de- 
gree than in the adiabatic regime. The behavior of the 
steady-state entanglement in the presence of mechanical 
damping is demonstrated in Fig. 3. We see that in this 
case the optimal entanglement does not occur in the adia- 
batical regime. With increasing thermal phonon number 
nth stronger coupling strengths gj are needed to achieve 
the maximum entanglement. However, the robustness of 
the generated entanglement is obvious, as it can still be 
maintained for a relatively high mean thermal phonon 
number n t h = 100. Reaching the quantum ground-state 
of the vibrational modes is therefore not a prerequisite of 
the present scheme, which reduces experimental difficul- 
ties considerably. 



III. MULTIPARTITE MECHANICAL 
ENTANGLEMENT 

In this section, we generalize the previous two-cavity 
model to a system of N mechanical oscillators in coupled 
cavities, as illustrated in Fig. 4 and proceed to discuss 
the generation of multipartite mechanical entanglement. 
Staying consistent with the previous model, we introduce 
the convention that the cavity modes a2 n -i and bi n are 
pumped by lasers which are blue-detuned from their res- 
onances, while ai n and 62™- 1 are driven by red detuned 
pumps. Therefore, the red and blue detuned modes are 
a2n-i and 62^-1 respectively for odd cavities, while they 
are bi n and a,2 n for even cavities, see Fig. l.(b). As- 
suming identical mechanical frequencies for all oscillators 



s=2 



cj = - J m Cj - ig a a* - ig b bj + y/2j^cf(t), 
where the symbols are 



(18a) 
(18b) 



(3a,3&,C?> C 7> a ?> 6 ?) = 



(51 , 92 , c\ , Cj , at , bj) for j odd 



(52,3l)Cj,c],a 5 -,6t) for j 



even 



with matched optomechanical couplings g a2n -i — 9b 2n — 
31 and 3f )2 „_ 1 = g a2n — 92, local vacuum noise operators 
z s (t), optical and mechanical loss-rates designated by k z 
and 7 m respectively and finally the coupling efficiencies 
rj z . Before we turn to the numerical solutions, let us 
first consider the situation that k z 3> {gj, 7m"-th}, which 
allows us to adiabatically eliminate the cavity modes. 
For the perfect intercavity couplings r\ z — 1 and identical 
cavity loss rates n z — k, the equations of motion for the 
odd and even mechanical oscillators are 



n-l 



C2n-l(t) 



-(7m + 7m)c2n-i(*) - 2jm ^2 C 2s -l(t) 

(19a) 



s=l 



v^^-iW+sm 



n-l 



C2n(t) = -(7™ + lm)C2n{t) - ^ C 2s (t) 



+ ^/^4 n n(^) + ^?(^)■ 



(19b) 



From the above equations we see that the odd and even 
mechanical oscillators are coupled to the noise operators 
C]"(t) and c!, n (i), respectively. Therefore, the entangle- 
ment may be established between any odd and even me- 
chanical oscillators with the nonclassical correlations be- 
tween the noises c^(t) and c™{t) given in Eq. (fT2l . How- 
ever, between oscillators with same parity, quantum en- 
tanglement cannot be established. This because that the 
source of entanglement in this scheme results from the 
coupling of the red sideband output into the blue side- 
band input and vice versa. For two even or odd oscil- 
lators, the cavity modes coupled to these two oscillators 



6 



0.6 
Ei: 



0.3 



-K-l 
K=5 
-K=10 



(a) 



0.07 



5 10 15 

gs/gi 





(c) 







E23 


m °2 




0.00 








A.1-23 






-0.04 




A.2-13 










XJ-12 






-0.08 






-0.12 




(b) 








-0.16 


\/ (d) . 



5 10 15 

g2/gi 



FIG. 5: The steady-state reduced bipartite entanglement E\2 
between the mechanical modes ci and C2 in (a) and the re- 
duced bipartite entanglement E23 between the modes 02, and 
C3 in (b) for a three-mode mechanical system. The tripartite 
entanglement structure of the mechanical oscillators is shown 
in (c) where the solid lines represent the reduced bipartite en- 
tanglement. The plot in (d) depicts the negative eigenvalues 
of the partially transposed correlation matrix of the three- 
mode mechanical system with respect to any one mechanical 
mode, which demonstrates that the genuine tripartite me- 
chanical entanglement can be achieved. The unit of k is 10 5 
Hz, the coupling strength gi/2n = 0.01 x 10 Hz, the mechan- 
ical rate 7 m /27r = 10 Hz, the mean thermal phonon number 
nth = 0, and the intercavity coupling efficiency rj — 1.0. 



have the same detunings from the pump lasers, which 
leads them not to being entangled but to the mode cou- 
pling through an incoherent exchange interaction with 
rate — 2"/ m . These results are verified in the following via 
numerical solution of Eqs. dTSl) . 

For AT = 3, we plot in Fig. 5 (a) and (b) the bipartite en- 
tanglement E\2 between the mechanical modes c\ and C2 
and the bipartite entanglement E23 of the modes C2 and 
C3, respectively. We see that the entanglement E\2 > E23 
for the same parameters. As predicted above, bipartite 
entanglement between the mechanical modes ci and C3 
is absent. Nevertheless, as demonstrated in Fig. 5 (d), 
full inseparable (genuine) tripartite entanglement can be 
established among the three remote mechanical oscilla- 
tors. Fig. 5 (d) depicts the negative eigenvalues A;_ mn 
of the partially transposed three-mode correlation ma- 
trix with respect to the l-th mode. The appearance of 
negative eigenvalue confirms bipartite entanglement be- 
tween the transposed mode I and the subsystem of the re- 
maining modes m and n, and fully inseparable (genuine) 
multipartite entanglement is demonstrated in the regime 
where the neg ative eigenvalues simultaneously exist for 
I = 1,2,3 [39(. Also, Fig. 5 (d) shows that bipartite en- 
tanglement between the mode C2 and the remaining two 
modes c\ and C3 is largest, since it is the only mode which 
is simultaneously entangled to the two other subsystems 
ci and C3. Finally, the entanglement between the mode 
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FIG. 6: (a) Steady-state bipartite entanglement between the 
mechanical modes C3 and C4 for a four-mode mechanical 
system, (b) Quadripartite square graph-state entanglement 
among the four remote mechanical oscillators, (c) Depen- 
dence of the reduced bipartite entanglement on the intercav- 
ity coupling efficiency r\. Other parameters used are the same 
as in Fig. 4. 



C3 and the subsystem including c\ and C2 is smallest, 
since the bipartite entanglement satisfies E23 < £-12 and 
E13 = 0. Therefore see that the negativities will satisfy 

A3— 12 > Al-23 > A2-13. 

When extending the above three-mode mechanical sys- 
tem to the four-mode case, i.e. N = 4, it is not difficult 
to see from Eq. (fl~8|) that the reduced bipartite entangle- 
ment E\2 and E23 are not affected due to the unidirec- 
tional cavity coupling. Therefore, the entanglement E12 
and E23 are the the same as in the N = 3 case plotted in 
Fig. 5. Furthermore, it can be inferred from Eg. (1191) that 
the bipartite entanglements will satisfy En — E23 in the 
bad-cavity limit. In Fig. 6 (a), we plot the bipartite en- 
tanglement E34 between the mechanical modes C3 and C4, 
and it is obvious that it exhibits similar behavior as the 
entanglement E12 between modes c\ and C2 ( see Fig. 5). 
We therefore see that quadripartite square graph-state 
entanglement among four remote mechanical oscillators 
can be achieved via cascaded cavity couplings. This kind 
of multipartite entanglement is useful in the field of long- 
distance quantum communication. The effects of imper- 
fect cavity couplings are illustrated in Fig. 6 (c). We see 
that for the coupling efficiency as low as rj = 0.5, the 
genuine quadripartite entanglement of four distant me- 
chanical oscillators can still be achieved. 



IV. CONCLUSION 

In conclusion, we propose a scheme to generate steady- 
state entanglement of remote mechanical oscillators in 
unidircctionally coupled cavities in the cascaded way. We 
note here that while the present model assumes the mem- 
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branes as mechanical oscillators, the role of mechani- 
cal elements can also be played by momentum modes 
of clouds of ultracold atoms. By choosing the detuning 
of the pump lasers, in each cavity the mechanical oscil- 
lator is coupled to the two cavity modes via parametric 
and beam-splitter-like interactions. The output quantum 
fluctuating field of the first cavity is subsequently driv- 
ing the second cavity with reversed detunings. For the 
case of two mechanical oscillators in cascaded cavities, 
the cavity dissipation can pull the two mechanical oscil- 
lators into a stationary two-mode squeezed vacuum state 
for negligible mechanical damping. The two-mode me- 
chanical entanglement depends on the relative strength 
of the pump lasers and is robust to thermal fluctuations. 
For multiple mechanical oscillators in multiple cascaded 
cavities, it is found that the steady-state bipartite entan- 
glement can be established between the odd and even os- 
cillators, whereas odd and even oscillators do not become 
entangled. We show that using this scheme the genuine 



multipartite entanglement can be achieved among remote 
mechanical oscillators by cavity dissipation. This kind of 
remote multipartite macroscopic entanglement is a use- 
ful resource in the construction of long-distance quantum 
communication networks. 
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